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Preface 


The  ^ '  o-dimensional  numerical  modeling  of  the  flow  conditions  in  the 
proposed  San  Timoteo  Channel  was  performed  for  the  U.S.  Army  Engineer 
District,  Los  Angeles  (SPL).  The  study  was  authorized  by  the  U.S.  Army 
Engineer  Division,  South  Pacific,  on  26  August  1993.  SPL  personnel 
involved  in  this  study  were  Mr.  Algis  Bliudzius  and  Ms.  Wendy  Gist  under 
the  supervision  of  Mr.  Brian  Tracy,  Chief  of  the  Hydraulics  Section. 

The  study  was  conducted  in  the  Hydraulics  Laboratory  of  the  U.S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  during  the  period  December 
1993  to  April  1994  uitder  the  direction  of  Messrs.  F.  A.  Herrmaiui,  Jr., 
Director  of  the  Hydraulics  Laboratory;  R.  A.  Sager,  Assistant  Director  of  the 
Hydraulics  Laboratory;  and  G.  A.  Pickering,  Chief  of  the  Hydraulic  Struc¬ 
tures  Division.,  Hydraulics  Laboratory. 

This  work  was  performed  by  Messrs.  Richard  L.  Stockstill  and  Mikel  W. 
Ott  under  the  super/ision  of  Mr.  John  F.  George,  Chief,  Locks  and  Conduits 
Branch,  Hydraulic  Structures  Division.  Technical  assistance  in  the  form  of 
consultation  and  peer  review  was  provided  by  Dr.  R.  C.  Berger,  Estuaries 
Division,  Hydraulics  Laboratory.  The  report  was  prepared  by  Mr.  Stockstill. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard.  EN 


The  coHleiUi  of  this  report  are  not  to  he  used  far  advertismg,  pubUcaliori, 
or  promolioual  purposes.  Chatiott  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows; 


MuMply 

By 

To  Obtain 

cubic  feet 

0.02831686 

cubic  meters 

degrees  (ertgle) 

0.01746329 

radians 

feat 


0.3048 


meters 


1  Introduction 


Background 

The  San  Timoteo  Creek  is  a  tributary  of  the  Santa  Ana  River  and  drains 
portions  of  San  Bernardino  and  Riverside  Counties  (Figure  1).  The  existing 
creek  has  the  capacity  to  protect  the  surrounding  community  from  approxi¬ 
mately  a  20-year-frequency  flood  (U.S.  Army  Engineer  District,  Los  Angeles, 
1993).  The  recommended  plan  would  provide  a  100-year  level  of  protection 
(19,000  cfs).' 


Figure  1 .  Vicinity  and  lojation  maps 

The  proposed  design  within  the  reach  studied  includes  a  sediment  basin 
(the  last  of  a  series  of  basins  for  the  project),  a  concrete  weir  followed  by  a 
chute  having  converging  sidewalls,  a  compound  horizontal  curve  consisting  of 
spirals  between  a  circular  curve  and  the  upstream  and  downstream  tangents 


'  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  found  on 
page  V. 
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with  a  banked  invert,  and  a  bridge  pier  associated  with  the  San  Timoteo 
Canyon  Road.  The  modeled  reach  extended  from  basin  sta  5+00  to  channel 
sta  280+00,  resulting  in  a  total  simulated  channel  length  of  2,396.94  ft. 
Geometrical  details  of  the  proposed  channel  design  are  shown  in  the  plan  view 
provided  in  Plate  1 , 


Purpose 

This  investigation  was  initiated  because  there  was  concern  as  to  the 
adequacy  of  a  one-dimensional  analysis  of  the  flow  conditions  within  the  chan¬ 
nel  chute.  There  was  a  question  as  to  whether  computing  cross-sectional  aver¬ 
aged  flow  variables  provided  a  sufficiently  accurate  estimate  of  flow  depths 
within  the  geometrically  complex  chute.  It  was  known  that  the  supercritical 
flow  in  the  chute  would  produce  oblique  standing  waves.  These  waves  do  not 
present  a  problem  as  long  as  the  channel  walls  are  of  sufficient  height.  How¬ 
ever,  the  flow  conditions  in  the  upper  end  of  the  channel  are  further  compli¬ 
cated  by  the  fact  that  the  chute  is  followed  by  a  channel  curve,  followed  by 
the  San  Timoteo  Canyon  Road  bridge.  A  two-dimensional  analysis  was 
deemed  necessary  to  evaluate  the  chute’s  influence  on  the  flow  conditions  in 
the  curve  and  the  curve’s  impact  on  the  flow  conditions  at  the  bridge. 


Approach 

The  two-dimensional,  depth-averaged  flow  model,  HIVEL2D,  was  used  to 
simulate  the  flows  in  the  high-velocity  channel.  This  model  was  chosen 
because  of  its  ability  to  simulate  supercritical  flow  and  capture  shocks  such  as 
oblique  standing  waves. 
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2  Model  Description 


Governing  Equations 

HIVEL2D  (Stockstill  and  Berger  1994)  solves  the  depth-averaged  unsteady 
shallow-water  equations  implicitly  using  finite  elements.  The  shallow-water 
equations  are  the  result  of  vertical  integration  of  the  equations  of  mass  and 
momentum  conservation  for  incompressible  flow  assuming  hydrostatic  pres¬ 
sure  distribution.  The  dependent  variables  of  the  two-dimensional  fluid 
motion  are  defined  by  the  flow  depth  and  the  x-direction  and  y-direction 
components  of  unit  discharge.  These  variables  are  functions  of  the  spatial 
Cartesian  coordinates  x  and  y  and  time.  A  brief  description  of  the  finite  ele¬ 
ment  model  is  provided  in  Appendix  A. 


Assumptions  and  Limitations 

An  important  (and  standard)  assumption  made  in  the  derivation  of  the 
governing  equations  is  that  vertical  accelerations  are  negligible  when  com¬ 
pared  to  the  horizontal  accelerations  and  the  acceleration  due  to  gravity.  This 
assumes  that  the  pressure  distribution  is  hydrostatic.  The  hydrostatic  pressure 
assumption  is  generally  applicable  to  high-velocity  channels;  however,  certain 
limited  areas  of  the  flow  field  are  decidedly  nonhydrostatic.  One  must  resort 
to  engineering  judgment  to  assess  the  ramifications  of  this  difference,  noting 
that  hydrostatic  models  overestimate  the  speed  of  short  wavelengths. 

Another  significant  assumption  is  that  the  invert  slope  is  geometrically  mild 
(less  than  about  0.02),  though  it  may  be  hydraulically  steep,  producing  super¬ 
critical  flow.  This  is  a  reasonable  assumption  for  most  high-velocity  channels. 
However,  the  flow  speed  in  long  channel  reaches  with  a  favorable  slope 
greater  than  about  O.OS  will  tend  to  be  overestimated  and  the  flow  depth  will 
tend  to  be  underestimated,  whereas  adverse  steep  slopes  will  tend  to  over¬ 
estimate  the  flow  depth  and  underestimate  the  flow  speed.  A  general  discus¬ 
sion  of  the  assumptions  associated  with  the  shallow-water  equations  is 
presented  in  Berger  and  Stockstill  (1994), 
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3  Model  Applications 


Application  of  the  model  to  the  proposed  channel  design  requires  the  con¬ 
struction  of  a  numerical  model  computational  mesh  to  represent  the  design's 
geometry  and  the  specification  of  boundary  conditions  and  model  parameters. 


Geometry  and  Computational  Mesh 

The  reach  of  the  San  Timoteo  Creek  channel  modeled  in  this  study  was 
from  basin  sta  5+00  to  0+20  and  from  channel  sta  299+ 16.94  to  280+00,  a 
total  length  of  2,396.94  ft.  The  subcritical  flow  in  the  sediment  basin  at  the 
upper  end  is  controlled  by  the  concrete  weir,  which  begins  the  concrete-lined 
high-velocity  channel.  The  upstream  face  of  the  weir  is  a  IV  on  2H  slope. 

The  weir  crest,  which  is  immediately  followed  by  a  940-ft  converging  sidewall 
chute,  is  20  ft  long  in  the  direction  of  flow.  The  chute  was  designed  to  transi¬ 
tion  from  the  338-ft-wide  sediment  basin  to  the  86-ft-wide  high-velocity  chan¬ 
nel.  A  channel  bend  is  located  83.35  ft  downstream  of  the  chute  end.  The 
bend  was  designed  with  spiral  curves  and  a  banked  invert.  The  San  Timoteo 
Canyon  Road  crosses  the  channel  downstream  of  the  bend.  One  pier  for  this 
bridge  is  located  within  the  channel,  near  the  channel  center  line.  Plate  1 
shows  a  plan  view  of  the  modeled  reach. 

The  width  of  the  flow  domain  within  the  sediment  basin  is  dependent  on 
the  water-surface  elevation  in  this  vicinity  because  the  channel  sidewalls  are 
sloping.  Therefore,  tlie  numerical  model  computational  mesh  of  the  sediment 
basin  was  different  for  the  two  discharges  simulated.  Initial  numerical  tests 
were  conducted  to  determine  an  adequate  top  bank  elevation  to  use  for  the 
domain  limit.  By  trial,  a  top  bonk  of  elevation  1316.0'  was  found  appropri- 
-ite  for  modeling  a  discharge  of  19,000  cfs.  Representing  the  basin  flow  using 
a  fixed  grid  was  deemed  adequate  because  the  flow  in  this  region  is  subcriti¬ 
cal.  Fixed  grid  boundary  conditions  are  sufficient  for  simulation  of  steady 
subcritical  flow  in  channels  having  sloping  sidewall  boundaries. 


'  Ail  elevations  (el)  cited  herein  are  in  feet  referenced  to  the  National  Geodetic  Vertical 
Datum  (NGVD). 
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The  San  Timoteo  Canyon  Road  bridge  pier  presently  has  a  semicircular 
nose  and  tail.  The  pier  was  modeled  with  a  triangular  nose  and  tail.  Signifi¬ 
cantly  more  grid  refinement  would  be  needed  to  simulate  a  semicircular  nose 
and  tail  and  was  not  deemed  necessary  for  this  study.  Proposed  channel 
modifications  include  the  addition  of  a  debris  nose  on  this  pier.  A  sloping 
debris  nose  cannot  be  accurately  modeled  with  a  depth-averaged  model;  there¬ 
fore,  only  the  existing  pier  configuration  was  simulated  in  this  study. 

The  numerical  model  computational  mesh  is  shown  in  Plate  2.  The  mesh 
consisted  of  945  nodes  and  918  elements.  Grid  resolution  was  more  refined 
in  the  vicinity  of  sidewall  boundary  alignment  changes  and  at  all  grade  breaks. 


Boundary  Conditions 


This  study  included  the  simulation  of  two  discharges,  the  design  discharge 
of  19,000  cfs  and  a  lesser  discharge  of  12,000  cfs.  The  inflow  boundary 
located  in  the  sediment  basin  was  subcriiical  and  therefore  required  two 
boundary  conditions.  The  x-  and  y-direction  components  of  unit  discharge 
were  specified  for  the  inflow  boundary  conditions.  These  unit  discharges 
were  computed  assuming  a  constant  velocity  over  the  inflow  boundary  cross 
section.  A  water-surface  elevation  at  the  inflow  boundary  was  assumed,  then 
the  inflow  velocity  was  computed  as  the  discharge  divided  by  the  flow  area. 
The  unit  discharges  were  calculated  as  the  product  of  the  flow  depth  at  each 
inflow  node  times  the  computed  velocity.  These  unit  discharges  were  further 
resolved  into  x-  and  y-components  normal  to  the  inflow  boundary.  The  out¬ 
flow  boundary  was  supercritical  and  therefore  required  no  boundary 
conditions. 


Model  Parameters 

Model  and  flow  parameters  used  in  the  initial  simulation  are  provided  in 
the  following  tabulation,  where  a  is  the  temporal  diffeiencing  weight  (a  =  1 
is  first-order  and  a  =  1.5  is  second-order  backward  differencing);  0^^  and  0^,, 
are  the  Petrov-Galerkin  parameters  for  smooth  flow  and  for  shocks,  respec¬ 
tively;  and  Qh  are  the  coefficients  used  in  determination  of  eddy 
viscosity;  and  g  is  acceleration  due  to  gravity.  Descriptions  of  these 
conditions  are  presented  in  Appendix  A, 


Modal  Condition 

Value 

a 

1.0 

^SM'  ^SH 

0.1,  0.5 

0.1,  0.6 

g 

32.2  ft/sec’ 
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Because  a  steady-state  solution  was  sought,  the  temporal  difference  order  was 
irrelevant.  Therefore,  a  first-order  lemporal  differencing  was  used  in  all 
simulation  runs. 
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4  Results 


As  stated  previously,  this  study  included  the  simulation  of  two  discharges, 
the  design  discharge  of  19,000  cfs  and  a  lesser  discharge  of  12,000  cfs.  This 
12,000-cfs  discharge  was  selected  for  testing  because  a  one-dimensional  analy¬ 
sis  conducted  by  the  U.S.  Army  Engineer  District,  Los  Angeles,  found  that 
the  12,000-cfs  discharge  actually  resulted  in  larger  Froude  numbers  in  the 
channel  chute.  It  was  felt  that  these  larger  Froude  numbers  might  be  accom¬ 
panied  by  higher  standing  waves  at  a  particular  station  than  those  produced  by 
the  design  discharge  (19,000  cfs). 

The  prototype  flow  geometry  and  boundary  roughness  for  the  concrete 
channel  reach  will  most  likely  result  in  a  Manning’s  n  value  of  approximately 
0.014.  The  concrete  reach  of  channel  extends  from  the  toe  of  the  upstream 
face  of  the  weir  (basin  sta  04-44  at  center  line)  downstream  to  the  end  of  tlie 
modeled  reach  (sta  2804-00).  However,  an  envelope  of  possible  hydraulic 
conditions  was  established  by  also  performing  simulations  of  the  flow  condi¬ 
tions  resulting  from  a  very  smooth  boundary  by  using  an  n  value  of  0.012. 

All  simulations  used  an  n  of  0.030  for  the  sediment  basin,  which  has  rip- 
rapped  side  slopes  and  a  natural  bed. 

The  sensitivity  of  simulation  results  to  the  choice  of  model  parameters  and 
mesh  refinement  was  determined.  This  sensitivity  was  established  by  chang¬ 
ing  one  condition  per  run  and  then  comparing  the  results  with  previous  runs  to 
assess  the  condition’s  impact  on  simulation  results. 


Design  Discharge,  19,000  cfs 

Initial  tests  (n  =  0.014) 

Water-surface  profiles  resulting  from  a  discharge  of  19,(XX)  cfs  and  an  n 
value  of  0.014  are  presented  in  Plate  3  and  in  tabular  form  in  Table  1 .  The 
profiles  illustrate  the  variation  in  the  depth  of  flow  at  the  sidewalls  resulting 
from  supercritical  flow  in  the  channel  chute.  A  significant  setup  is  noted  in 
the  channel  bend.  It  is  interesting  to  note  the  rise  in  the  water-surface  eleva¬ 
tion  in  the  basin  as  the  flow  approaching  the  structure  decelerates.  Of  course 
the  total  energy  grade  line  decreases  in  the  direction  of  flow,  but  due  to  the 
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gradual  expansion  of  the  basin,  the  velocity  head  reduction  is  greater  than  the 
loss  of  total  head.  It  is  clearly  shown  in  the  profile  that  sufficient  clearance  is 
provided  through  the  San  Timoieo  Canyon  Road  bridge.  The  significant 
runup  on  the  upstream  face  of  the  pier  is  more  of  a  model  artifact  than  what 
would  actually  occur.  Vertical  accelerations  immediately  upstream  of  the  pier 
nose  are  of  such  magnitude  that  the  hydrostatic  assumption  is  not  valid. 

A  hydrostatic  model  provides  a  poor  estimate  of  this  sudden  runup.  To 
illustrate  this  point,  a  first  approximation  of  the  magnitude  of  these  accelera¬ 
tions  results  in  the  conclusion  that  this  pier  nose  runup  is  a  violation  of  the 
hydrostatic  assumption  and  that  the  runup  is  overestimated.  A  first  approxi¬ 
mation  of  the  magnitude  of  these  accelerations  is  obtained  by  applying  the 
kinematic  boundary  condition  to  the  water  surface: 


w 


i 


■57 


^  *  ivii^ 

dt  '  ‘  as 


(1) 


where 

w,  =  vertical  component  of  velocity  at  the  water  surface 
h  =  flow  depth 
t  =  time 

s  =  stream  direction 
and 


V|  =  (M^  +  V^) 


1/2 


where  «  and  v  are  the  x-  ?,nd  y  direction  components  of  velocity,  respectively, 
The  vertical  acceleration  a,  is: 


'  dt  dt  '  '  ds  dz 


(2) 


Neglecting  the  last  term  and  applying  to  steady  flow, 


I  I/I 

'''III 


(3) 


and 
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(4) 


'  d* 


I VI 


ds 


VI  ^ 

as 


Using  difference  approximations  of  the  results  of  a  simulation  of 
19,000  cfs  and  an  n  value  of  0.014,  the  ratio  of  a^/g  at  the  pier  nose  was 
calculated  to  be  9.7.  This  definitely  violates  the  hydrostatic  assumption,  i.c., 
that  the  vertical  accelerations  are  negligible  when  compared  to  gravity  (Appen¬ 
dix  A)  Pressures  at  the  pier  nose  in  the  actual  system  (prototype)  will  be 
greater  than  hydrostatic,  which  will  act  to  suppress  the  water-surface  rise. 
Determination  of  what  the  water-surface  elevation  will  be  in  the  prototype  is 
difficult.  However,  it  will  be  significantly  less  than  that  calculated  assuming  a 
hydrostatic  pressure  distribution. 


To  illustrate  that  the  hydrostatic  assumption  is  appropriate  elsewhere  in  the 
flow  field,  the  vertical  acceleraiion  was  approximated  in  the  area  of  the 
oblique  standing  wave  initiated  at  the  pier  nose  (sta  282  +  10,  10  ft  right  of 
center  line).  The  vertical  acceleration  in  this  vicinity,  relative  to  gravity,  was 
calculated  to  be  0.4.  Although  this  is  not  insignificant,  it  shows  that  the 
hydrostatic  assumption  is  reasonable  even  in  areas  where  the  flow  is  extremely 
rough. 


Depth  Contours  for  a  discharge  of  19,000  cfs  and  an  n  value  of  0.014  are 
presented  in  Figure  2.  Depth  contours  within  the  chute  are  shown  in  Fig¬ 
ure  3.  The  contours  illustrate  the  oblique  standing  waves  resulting  from  the 
converging  sidewalls.  Vertically  exaggerated  water-surface  mesh  plots 
(Plates  4  and  S)  .show  the  variation  in  depths  and  the  superelevation  in  the 
channel  bend.  Depth  contours  in  the  vicinity  of  the  bridge  are  provided  in 
Figure  4. 

At  this  point  in  the  study,  there  was  interest  in  the  effect  of  the  oscillations 
of  the  computed  depths  at  the  sidewalls  in  the  vicinity  of  the  crest  on  the 
results  dowiiStream,  particularly  if  these  oscillations  were  a  numerical  rather 
than  a  physical  phenomenon.  A  second  run  using  the  discharge  of  19,0(X)  cfs 
and  an  n  value  of  0.014  was  made  using  a  different  Petrov-Galerkin  weighting 
parameter  for  smooth  flow,  This  simulation  used  a  of  0  25  rather 
than  aPs^ofO.!,  which  was  used  in  the  initial  simulation.  The  0,^  is  the 
minimum  damping  parameter  used  throughout  the  flow  field. 

The  results  obtained  from  the  two  /3sm’s  are  compared  in  Plate  6.  These 
profiles  show  that  the  variation  of  0^  resulted  in  only  minor  changes  in  the 
computed  depths.  Moreover,  the  two  runs  converged  to  practically  the  same 
solution  150  ft  downstream  of  the  crest.  Therefore,  it  is  apparent  that  the 
solution  of  this  particular  flow  field  is  not  significantly  influenced  by  the 
choice  of  the  Petrov-Galerkin  weighting  parameter, 

The  model’s  simulation  of  the  flow  in  the  crest  vicinity  was  further 
analyzed  by  grid  refinement.  The  density  of  the  original  mesh  was  doubled  in 
both  directions  on  and  adjacent  to  the  crest.  Details  of  the  original  and 
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Figure  2.  Depth  contours  for  the  design  discharge  (n  =  0.014) 
refined  grids  in  this  area  are  provided  in  Plate  7. 

Water-surface  profiles  for  the  design  discharge  and  an  n  value  of  0.014 
computed  using  the  refined  grid  are  shown  in  Plate  8  and  Table  2.  The 
refined  solution  also  oscillated  at  the  crest.  The  results  obtained  from  the 
original  and  fine  grids  =  0.1  for  both  tests)  are  compared  in  Plate  9. 
Again  it  is  apparent  that  a  node-to-node  oscillation  occurred  upstream  and 
downstream  of  the  point  at  which  the  flow  was  accelerating  from  subcritical  to 
supercritical  regardless  of  grid  refinement.  The  oscillations  were  pronounced 
only  at  the  sidewalls  where  the  flow  passed  around  the  wing  walls  and 
intercepted  the  converging  sidewalls.  Therefore,  it  is  difficult  to  determine 
whether  these  results  are  a  numerical  artifact  or  physical.  However,  the 
solutions  converged  at  a  distance  of  150  ft  downstream  of  the  crest. 

Therefore,  using  the  original  grid  and  a  of  0.1  are  adequate  for  design 
purposes  and  were  used  throughout  the  remainder  of  the  study. 
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Figure  4.  Depth  contours  in  vicinity  of  San  Timoteo  Canyon  Road  bridge  pier  (discharge 
19,000  cfs  and  n  =  0.014) 


elevations  eitmputcd  using  an  n  value  of  0.012  were  less  than  or  equal  to 
those  computed  using  an  n  value  of  0.014. 


12,000-cfs  Discharge 

A  new  numerical  model  computational  mesh  was  constructed  for  the  simu¬ 
lation  of  a  discharge  of  12,(KK)  cfs.  This  new  mesh  was  needed  because  the 
domain  width  in  the  sediment  basin  was  reduced  from  that  resulting  from  a 
discharge  of  10,(K)0  cfs.  The  top  bank  elevation  was  selected  by  first  deter¬ 
mining  the  expected  water-.surface  elevation  in  the  basin  for  a  discharge  of 
l2.(K)0  cfs.  Given  the  water-surface  elevation  simulated  for  a  discharge  of 
19,000  cfs,  a  discharge  coefficient  for  the  "broad-crested  weir"  was  cal¬ 
culated,  This  discharge  coefficient  was  then  used  to  compute  the  basin  water- 
surface  elevation  for  a  di.scharge  of  12,000  cfs.  It  was  then  determined  that  a 
top  bank  elevation  of  1314.0  was  appropriate  to  use  to  construct  the  grid  for 
the  lesser  discharge.  The  numerical  model  computational  mesh  for  a  dis¬ 
charge  of  12,000  cfs  differed  from  the  original  mesh  only  in  the  vicinity  of 
the  side  boundaries  within  the  basin. 


Simu)ati'’n  results  with  a  discharge  of  12,000  cfs  for  n  values  of  0  014  and 
0.012  are  pre.sented  as  water-surface  profiles  in  Plates  1 1  and  12  and  Tables  4 
and  5,  respectively.  As  expected,  the  depths  are  less  than  those  produced  by 
the  design  discharge  and  no  adverse  standing  waves  were  noted.  The 
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maximum  water-surface  elevations  resulting  from  a  discnarge  cf  12,000  cfs 
using  either  an  n  value  of  0.014  or  an  n  value  of  0.012  were  less  than  those 
computed  for  the  design  discharge  with  an  n  value  of  0.014. 
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5  Summary  and  Conclusions 


Two-dimensional  (depth-averaged)  simulations  of  the  flow  conditions  in  the 
San  Timoteo  Channel  were  conducted  to  assess  the  adequacy  of  the  proposed 
design.  Particular  attention  was  given  to  the  variation  in  flow  depths  within 
the  entrance  chute.  The  flow  at  the  crest  accelerated  from  subcritical  to 
supercritical  flow  and  remained  supercritical  through  the  chute  and  channel 
thereafter.  Oblique  standing  waves  were  produced  when  the  supercritical  flow 
intercepted  the  chute's  laterally  converging  sidewalls.  However,  these  stand¬ 
ing  waves  were  expected  and  given  that  the  sidewalls  will  be  designed  of  suf¬ 
ficient  height,  the  flow  conditions  were  deemed  acceptable. 

The  prototype  flow  geometry  and  boundary  roughness  will  most  likely 
result  in  a  Marming’s  n  value  approximately  equal  to  0.014.  However,  flow 
conditions  resulting  from  a  Manning’s  n  value  of  0.012  were  also  simulated  so 
that  an  envelope  of  possible  hydraulic  conditions  could  be  established.  Along 
this  same  line  of  reasoning,  in  addition  to  the  design  discharge  (19,0(X)  cfs),  a 
smaller  discharge  (12,(K)0  cfs)  was  simulated.  For  a  given  discharge,  there 
was  no  significant  difference  in  the  water-surface  elevations  computed  using 
an  n  value  of  0.014  and  an  n  value  of  0.012.  The  smaller  n  value  did,  how¬ 
ever,  result  in  the  wave  crests  being  located  further  downstream  due  to  the 
increase  in  velocity.  No  adverse  flow  conditions  were  noted  for  a  discharge 
of  12,000  cfs. 

Simulation  results  indicate  that  the  proposed  San  Timoteo  Channel  design 
and  in  particular,  the  San  Timoteo  Canyon  Road  bridge,  will  convey  the 
design  discharge  of  19,000  cfs  in  an  acceptable  manner.  The  prototype  will 
experience  significantly  less  runup  on  the  pier  nose  than  the  simulation  results 
indicated.  Sidewall  heights  should  be  designed  based  on  the  simulation  results 
realizing  the  uncertainties  associate  with  the  analyses. 
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298  +  86.94 

1316,00 

298  +  96.94 

1316.91 

296  +  96.94 

1312,09 

298  +  81,94 

1315.72 

298  +•  91.94 

1316.36 

296  4-  66.94 

1311.64 

298  +  76,94 

1316.65 

298  +  86.94 

1316.19 

296  4-  36,94 

1311.06 

298  +  66.94 

1315.30 

298  +  81.94 

1315.96 

296  4-  06.94 

1310.52 

298  +  46,94 

1314.73 

298  +  76.94 

1315,78 

295  4-  74.94 

1309,96 

298  +  22.94 

1314.21 

298  +  66.94 

1316.61 

296  4-  46.61 

1309.74 

297  +  90.94 

1313.65 

298  +  46.94 

1314.93 

296  +  16.28 

1308.87 

297  4  68.94 

1312,94 

298  +  22.94 

1314.38 

294  +  86.94 

1308.44 

297  4  26.94 

1312,36 

297  +  90.94 

1313.69 

294  +  46.94 

1308.54 

296  4  96.94 

1311.68 

297  +  58.94 

1313.06 

294  4-  06.94 

1308,15 

296  4  66.94 

1311 .23 

297  +  26.94 

1312,47 

293  4-  64.94 

1307.28 

296  4  36.94 

1310.29 

296  f  96.94 

1311.81 

293  4-  22.94 

1305.23 

296  4  06,94 

1309.67 

296  +  66.94 

1311.27 

292  +  80.94 

1302.88 

295  4  74.94 

1310.33 

PU-atiAMunBiB— 

292  +  38.94 

1302.03 

295  4  46.61 

1310.19 

296  +  06.94 

1309.89 

291  +  96.94 

1302.10 

295  4  16.28 

1310.00 

295  +  74.94 

1310.23 

291  +  64.94 

1301 .71 

294  4  86.94 

1309.63 

295  +  45.61 

1310.11 

291  4-  12.94 

1299.83 

294  4  46.94 

1308.39 

295  +  16.28 

1309.88 

290  4-  70.94 

1297.87 

294  4  26,94 

1308.08 

294  +  86.94 

1309.60 

290  4-  28.94 

1297.20 

294  4  06.94 

1306,88 

294  +  46.94 

1308,30 

290  +  07.94 

1297.13 

293  4  85.94 

1304.85 

(Continued) 


Note:  Sides  of  channel  are  referenced  to  looking  downstream. 


able  2  (Concluded) 


tation 


294  +  26.94 


294  -I-  06.94 


293  +  85.94 


293  +  64.94 


293  +  22.94 


292  +  80.94 


292  +  38.94 


291  +  96.94 


291  +  54.94 


291  +  12.94 


290  +  70.94 


290  +  28.94 


290  +  07.94 


289  +  86.94 


289  +  66.29 


289  +  45.64 


289  +  03.59 


288  +  70.26 


288  +  36.92 


288  +  03.59 


287  +  61.97 


287  +  20.35 


286  +  78.74 


286  +  37.12 


285  +  9B.50 


285  +  53.88 


285  +  12.26 


284  +  70.64 


284  -t  29.03 


283  +  87.41 


283  +  46.79 


283  +  12.46 


282  +  79.12 


282  +  62.46 


282  +  45.79 


282  +  32.98 


282  +  20.12 


282  +  09.1  1 


281  +  98.1  1 


281  +  87.10 


281  +  76.09 


281  +  59.61 


281  +  43.10 


281  +  27.20 


280  +  95.40 


280  +  63.60 


280  +  31.80 


280  +  00.00 


Elavation 
Left  Side 


tation 


289  -I-  86.94 


289  +  66.29 


289  +  45.64 


289  +  03.59 


288  +  70.26 


288  +  36.92 


288  +  03.59 


287  +  61.97 


287  +  20.35 


286  +  78.74 


286  +  37.12 


285  +  95.50 


285  +  53.88 


285  +  12.26 


284  +  70.64 


284  +  29.03 


283  +  87.41 


283  +  45.79 


283  +  12.46 


282  +  79,12 


282  +  62.46 


282  +  45.79 


282  +  38.08 


282  +  30.38 


282  +  25.25 


282  +  20,12 


281  +  76.12 


281  +  75.50 


281  +  71.89 


281  +  65.46 


281  +  59.22 


281  +  50.18 


281  +  39.32 


281  +  23.85 


280  +  95.40 


280  +  63.60 


280  +  31.80 


280  +  00.00 


tation 


293  +  64.94 


293  +  22.94 


292  +  80,94 


292  +  38.94 


291  +  96.94 


291  +  54.94 


291  +  12.94 


290  +  70.94 


290  28.94 


290  +  07.94 


289  4-  86.94 


289  4-  66.29 


289  4-  45.64 


289  4-  03.59 


288  4-  70.26 


288  4-  36.92 


288  4-  03.59 


287  4-  31.97 


287  4-  20.35 


286  4-  78.74 


286  4-  37.12 


285  4-  95,50 


285  4-  53.88 


285  4-  12.26 


284  4-  70.64 


284  4-  29.03 


283  4-  87.41 


283  4-  45.79 


1308.01 


1306.76 


1304,71 


1305.23 


1305.40 


1304.67 


1303.36 


1301,38 


1299.88 


1299.67 


1299.05 


1297.80 


1297.24 


1296.15 


1294.92 


1294.74 


1295.14 


1294.07 


1292.74 


1291.14 


1290.64 


1290.67 


1290.25 


1289.32 


1288.38 


1287.73 


1287.30 


1286.94 


1286.55 


1286.14 


1285,63 


1286.01 


1286.51 


1286,72 


1286.84 


1286.43 


1286.23 


1286.13 


1285.95 


1286.03 


1286.35 


1285.92 


1286.97 


1286.75 


1285.23 


1285.10 


1285.90 


1285.36 


1297.09 


1296.88 


1296.38 


1294.99 


1293.92 


1293.39 


1 293.74 


1293.67 


1292.76 


1291.81 


1 291.32 


1291.06 


1290.67 


1290.06 


1289.44 


1288.89 


1288.43 


1287.89 


1287.49 


1286.98 


1286.64 


1286.47 


1286.44 


1286,17 


1286.47 


1295.00 


1282.46 


1281.62 


1233.48 


1284.66 


1285.33 


1285.84 


1286.07 


1285.85 


1285.85 


1286.11 


1285.78 


1285.18 


1305.26 


1305.32 


1304.58 


1303.28 


1301.27 


1299.90 


1299.83 


1299.13 


1297.73 


1297.52 


1296.11 


1294.81 


1294.55 


1295.26 


1296,03 


1295.63 


1295.48 


1294,72 


1294.29 


1293.94 


1293.58 


1293,14 


1292,74 


1292,31 


1231.74 


1291.05 


1290.38 


1289.83 


283  +  12.46 

1288.70 

282  +  79.12 

1287.05 

282  +  62.46 

1286.95 

282  +  45,79 

1286.76 

282  +  32.98 

1286.67 

282  +  20.12 

1286.53 

282  +  09.1  1 

1286,37 

281  +  98.11 

1286.59 

281  +  37.10 

1286.18 

281  +  76.09 

1286.33 

281  +  59,61 

1287,48 

281  +  43.10 

1287.34 

281  +  27.20 

280  +  95.40 

1284.96 

280  +  63.60 

1284.9 

280  +  31.80 

1285.31 

280  +  00.00 

1285.32 

able  3 

Vater-Surface  Elevations.  Discharge  19,000  ofs,  n  =  0.012 


5  +  00.00  11317.44 


4  +  42.82  1317.56 


3  +  83.61 


3  +  83.60 


3  +  43.97 


3  +  06.91 


2  +  70.02 


2  +  29.47 


1  +  87.11 


1  +  61,01 


1  +  39  81 


1  +  18.76 


0  +  97.88 


0  +  81.31 


0  +  68.04 


1317,64 


1317.50 


1317.67 


1317.79 


1317.81 


1317.87 


1317.77 


1317.76 


1317.84 


1317.82 


1317.84 


1317.83 


1317.79 


0  +  65.40  1317.86 


0  +  43.51 


0  +  31.61 


1317.63 


1318.02 


299  +  16.94 


299  +  06.94 


1317.91 


1316.48 


298  +  96.94  1317.39 


298  +  86.94 


298  +  76.94 


1315.89 


1316.00 


298  +  61,94  1316.29 


298  +  46.94  1314.94 


298  +  22.94  1314,29 


297  +  90.94  1313.64 


297  +  58.94  1312,98 


297  +  26.94 


296  +  96.94 


296  +  66.94 


296  +  36.94 


296  +  06.94 


295  +  74.94 


295  +  45.61 


295  +  16.28 


1312.42 


1311.69 


131 1.15 


1310.61 


1309.61 


1309.38 


1309.85 


1309.78 


294  +  86.94  11309.37 


2S4  +  46.94  1308.39 


1308.13 


294  +  06  94  11306.73 


293  +  85,94 


293  +  64.94 


1304.73 


1304.18 


1304.70 


6  +  00,00  1317.50 


4  +  40.00  1317.52 


3  +  80.00  1317.62 


3  +  20.00  1317.69 


2  +  72.00  1317.75 


2  +  24.00  1317.78 


1  +  88.00  1317,80 


1  +  64.00  1317.81 


1  +  40.00  1317.83 


5  +  00.00  1317.74 


1  +  16.00 


0  +  92.00 


0  +  44.00 


299  +  16.94 


298  +  96.94 


298  +  76.94 


29B  +  46.94 


298  +  22.94 


297  +  90.94 


294  +  86.94  1308.26 


1317,83 


1317.88 


0  +  68.00  1317.82 


1317.96 


1317.34 


1314.68 


1313.97 


1313.30 


1313.01 


1311.71 


297  +  58.94  1311.52 


297  +  26.94  1312.64 


296  +  96.94  1312.12 


296  +  66.94  1311.60 


296  +  36,94  1311.11 


296  +  06,94  1310.55 


296  +  74,94  1309.81 


296  +  45.61  1309,43 


295  +  16.28  1309.02 


1294  +  46.94  130’ 


1307.8/ 


|293  +  64,94  |1 307.68 


293  +  22.94 


292  +  80.94 


292  +  38.94 


291  +  96.94 


291  +  54.94 


231  +  12.94 


290  +  70.94 


290  +  28.94 


290  +  07.94 


1305.41 


1302.89 


:  301. 21 


1301.12 


1301.79 


1300.68 


1298.05 


1296.36 


1295.99 


|289  +  86.94  [1296.19 


289  66.29  |l  296.26 


|289  +  45.64  1296.82 


4  +  50.89 


3  +  91.28 


3  +  46.00 


2  +  87.60 


2  +  47.99 


2  +  1 1.60 


1  +  76.99 


1  +  65.55 


1  +  03.99 


0  +  79.99 


0  +  55.39 


0  +  42.45 


0  +  31.32 


299  +  16.94 


299  +  16.94 


299  +  06.94 


298  +  96.94 


298  +  86.94 


298  +  76.94 


298  +  61,94 


298  +  46.94 


298  +  22.94 


297  +  90.94 


297  +  58.94 


297  +  26.94 


296  +  96.94 


296  +  66.94 


1317.15 


1317.60 


1317.66 


1317.73 


1317.74 


1317.78 


1317.77 


1317.80 


1  +  27.99  1317,78 


296  +  06.94 


295  +  74.94 


295  +  45.61 


295  +  16.28 


294  +  86.94 


294  +  46.94 


294  +  26.94 


294  +  06.94 


293  +  85.94 


293  +  64.94 


293  +  22.94 


292  +  80.94 


292  +  38.94 


291  +  96.94 


1317.82 


1317.76 


1317.82 


1317,64 


1317.92 


1317.85 


1317.77 


1314.74 


1317.46 


1315.45 


1316.91 


1315.03 


1314.75 


1314.11 


1313.53 


1312.84 


1312.35 


1311.55 


131 1.04 


1310.62 


1309.40 


1309.26 


1309,97 


1309.86 


1309.46 


1308.63 


1308.23 


1306.84 


1304.79 


1304.28 


1304.77 


1304.49 


1303.63 


1301.63 


(Continued) 


ote:  Sides  of  channei  are  referenced  to  looking  downstream. 


able  3  (Concluded) 


II 

Elevation 

1 

tation 

Right  Side 

tation 


292  +  80,94 


292  +  38.94 


291  +  96.94 


291  +  54.94 


291  +  12.94 


290  +  70.94 


290  +  28.94 


290  +  07,94 


289  +  86.94 


289  +  66.29 


289  +  45.64 


289  +  03.59 


288  +  70.26 


288  +  36.92 


288  +  03.59 


287  +  61 .97 


287  +  20.36 


286  +  78.74 


286  +  37.12 


285  4-  96.50 


285  4-  53.88 


1304.62 


1303.72 


1301.68 


1299.45 


1298.81 


1298.93 


1298.00 


1297.39 


1296.35 


1294.28 


1293.95 


1294.81 


1293.77 


1292.81 


1291.05 


1290.06 


1289.79 


1289.65 


1289.01 


1288.02 


287  +  20,35 


286  4-  78.74 


286  4-  37.12 


285  4-  95.50 


285  4-  53.88 


285  4-  12.26 


284  4-  70.64 


284  4-  29.03 


283  4-  87.41 


283  4-  45.79 


283  +  12.46 


282  4-  79.12 


282  4-  62.46 


282  4-  45.79 


282  (•  38.08 


282  4-  30.38 


282  4-  25.26 


1295.65 


1293.85 


1292.55 


1292.74 


1293.40 


1292.98 


1291.83 


1290,83 


1290,44 


1290.33 


1289.96 


1289.28 


1288.52 


1287.97 


1287.51 


1287.20 


1286.73 


1286.39 


1286.17 


1286.12 


1285.85 


1285.15 


291  4-  64.94 


291  +  12.94 


290  +  70.94 


290  4-  28.94 


290  4-  07.94 


289  +  86.94 


289  4-  66.29 


289  4-  45.64 


289  4-  03.59 


288  4-  70.26 


288  4  36.92 


288  4-  03.59 


287  4-  61.97 


287  4-  20.36 


286  4-  78.74 


286  4-  37.12 


285  4-  95.50 


285  4-  53.88 


285  4-  12.26 


284  4-  70.64 


284  4-  29.03 


283  4-  87.41 


11299.34 


1298.91 


1299.06 


1298.00 


1297,35 


1296.28 


1294.50 


1293.76 


1295.04 


1295.79 


1295.60 


1295.55 


1294.57 


1293.96 


1293.86 


1293.70 


1293.18 


1292.57 


1292.06 


1291.60 


1291.08 


1290.44 


284  4-  70.64 

1286.33 

282  4-  20,12 

1294.72 

283  4-  45.79 

1289.81 

284  4-  29,03 

1286.02 

281  4-  76.12 

1282.46. 

233  4-  12.46 

1288.59 

283  4-  87.41 

1285.66 

281  4-  75.50 

1281.35 

282  4-  79.12 

1286.83 

283  +  46.79 

1284.94 

281  4  71  89 

1282.94 

282  4-  62.46 

(286.67 

283  4-  12.46 

1285.31 

281  4-  65.46 

1284,10 

282  4-  45,79 

1286.46 

282  4-  79.12 

1285.88 

281  4-  69.22 

1284,84 

282  4-  32.98 

1286.34 

282  4-  62.46 

1236.17 

281  4-  50.18 

1285.30 

282  4-  20.12 

1286.19 

281  4-  39.32 

1285.78 

282  +  09.1 1 

1286,02 

281  4-  23.85 

1286.72 

281  +  98.11 

1286.08 

282  4-  20.12 

280  +  95,40 

1285.43 

261  4-  87,10 

1286.10 

282  4-  09.1 1 

280  4-  63.60 

1285.57 

281  4-  76.09 

1285.47 

281  4-  98.1 1 

1285.72 

280  +  31.80 

1285.51 

281  4-  59.61 

1286.47 

281  4-  87.10 

1285.84 

280  4  00.00 

1285.03 

281  4-  43.10 

1287.08 

281  4-  43.10  1286.37 


281  4-  27.20  1286,53 


280  4-  95.40  128  1.99 


280  4-  63.60  1284.49 


280  4-  31 .80  1284.62 


280  4-  00.00  1284.86 


281  4-  27.20 


280  4-  95.40 


280  4-  63.60 


280  4-  31.80 


280  4-  00.00 


1286.41 


1284.80 


1284.37 


1284.59 


1284.71 


Table  4 

Water-Surface  Elevations,  Discharge  1 2,000  c 

fs. 

n 

=  0.014 

Station 

Elavation 

Laft  Sida 

Station 

Elavation 

Cantar 

Station 

5  +  00.00 

1315.82 

5  +  00.00 

1315.89 

5 

4 

00.00 

1315.99 

4  +  42.82 

1315.94 

4  +  40.00 

1315.90 

+ 

50.89 

1315,72 

3  +  83.61 

1315.93 

3  +  80.00 

1315.95 

3 

4 

91.28 

1315.95 

3  +  83.60 

1315.91 

3  +  20.00 

1315.99 

3 

4  46.00 

1315.97 

3  +  43.97 

1315.96 

2  +  72.00 

1316.02 

2 

+ 

87.60 

1316,01 

3  +  06.91 

1316.05 

2  +  24.00 

1316.04 

2 

4  47.99 

1316.01 

2  +  70.02 

1316.05 

1  +  88.00 

1316.05 

2 

+ 

11.60 

1316.04 

2  +  29.47 

1316.08 

1  +  64.00 

1316.05 

1 

+ 

75.99 

1316.03 

1  +  87.11 

1316.04 

1  4  40.00 

1316.06 

1 

+ 

55.56 

1316.05 

1  +  61.01 

1316.03 

1  4  16.00 

1316.06 

1 

+ 

27.99 

1316.03 

1  +  39.81 

1316.07 

0  4  92.00 

1316.09 

1 

4 

03.99 

1316.06 

1  +  18.76 

1316.05 

0  4  68,00 

1316.05 

0 

+ 

79.99 

1316.02 

0  +  97.88 

1316.07 

0  4  44.00 

1316.15 

0 

+ 

55.39 

1316.03 

0  +  81.31 

1316.05 

299  4  16.94 

1315.74 

0 

4 

42.45 

1315.94 

0  +  68.04 

1316.03 

298  4  96.94 

1313.61 

0 

4 

31.32 

1315.98 

0  +  65.40 

1316.07 

298  4  76.94 

1313.04 

299 

+ 

16.94 

1315.72 

0  +  43.51 

1315.94 

298  4  46.94 

1312,35 

299 

4 

16.94 

1315.78 

0  +  31.61 

1316.09 

298  4  22.94 

1312.00 

299 

4  06.94 

1314.50 

299  +  16.94 

1316.01 

297  4  90.94 

1311.32 

298 

+ 

96.94 

1316.44 

299  +  16.94 

1315.80 

297  4  68.94 

1310.41 

298 

4 

86.94 

1314.50 

299  +  06.94 

1314.54 

297  4  26.94 

1310.81 

298  4 

76.94 

1314,25 

298  +  96.94 

1315.54 

296  4  96.94 

1310.91 

298 

4  61.94 

1313.81 

298  +  86.94 

1314.58 

296  4  66.94 

1310,21 

298 

4  46,94 

1313.46 

298  +  76.94 

1314,40 

296  4  36.94 

1309.69 

298 

+ 

22.94 

1312.95 

298  +  61,94 

1313,91 

296  4  06.94 

1309,32 

297 

+ 

90.94 

1312.31 

298  +  46.94 

1313.57 

295  4  74.94 

1308,76 

297 

+ 

58.94 

1311.81 

298  +  22,94 

1313,03 

295  4  45.61 

1308.22 

297 

+ 

26.94 

1311.19 

297  +  90.94 

1312.40 

295  4  16.28 

1307.87 

296 

96.94 

1310.71 

297  +  58.94 

1311,85 

294  4  86.94 

1307.34 

296 

+ 

66,94 

1309.90 

297  +  26.94 

131 1.27 

294  4  46.94 

1306.1 1 

296 

+ 

36.94 

1309.60 

296  +  96,94 

1310.72 

294  4  06,94 

1305.15 

296 

■f 

06.94 

1308.89 

296  +  66.94 

1310.01 

293  4  64.94 

1305.51 

295 

4 

74.94 

1308.05 

296  +  36,94 

1309.61 

293  4  22.94 

1303.99 

295 

+ 

45.61 

1307,92 

296  +  06.94 

1308.91 

292  4  80.94 

1301,96 

295 

4 

16,28 

1308.12 

295  +  74.94 

1308.16 

292  4  38.94 

1300,15 

294 

4 

86.94 

1307.87 

295  +  45.61 

1308.01 

291  4  96.94 

1298.64 

294 

4  46,94 

1306.76 

295  +  16.28 

1308.06 

291  4  64,94 

1298.45 

294 

+ 

26.94 

1306.58 

294  +  86.94 

1307,83 

291  4  12,94 

1298,29 

294 

+ 

06.94 

1305.66 

294  +  46.94 

1306.74 

290  4  70,94 

1297.19 

293 

+ 

85.94 

1303.88 

294  +  26.94 

1306.52 

290  4  28.94 

1295.70 

293 

+ 

64.94 

1303.43 

294  +  06.94 

1305.58 

290  4  07.94 

1294.70 

293 

+ 

22.94 

1302.57 

293  +  85.94 

1303.78 

289  4  86.94 

1293.91 

292 

+ 

80.94 

1302.00 

29.3  +  64.94 

1303.47 

289  4  66,29 

1293.44 

292 

+ 

38.94 

1301.65 

293  +  22.94 

1302.49 

289  4  45.64 

1293.36 

291 

+ 

96.94 

1300.41 

(Continued) 


Note:  Sides  of  channel  are  referenced  to  looking  downstream. 


able  4  (Concluded) 


1301.88 

1301.73 

1300.45 

291  +  54.94 

1298.46 

291  +  12.94 

1297.18 

290  +  70.94 

1296.24 

290  +  28,94 

1295.63 

290  +  07,94 

1295.70 

289  +  86.94 

1294.42 

289  +  66.29 

1293.30 

289  +  45.64 

1292.70 

289  +  03.69 

1292.47 

288  +  70.26 

1291.57 

288  +  36.92 

1290.76 

288  +  03.69 

1289.58 

287  +  61.97 

1288.97 

287  +  20.35 

1288.79 

286  +  78.74 

1288.63 

286  +  37.12 

1288.09 

285  +  96.50 

1287.31 

285  +  53,88 

1286.60 

285  +  12.26 

1285.79 

284  +  70.64 

1285.22 

284  +  29.03 

1284.81 

283  +  87.41 

1284.53 

283  +  45,79 

1284,21 

1284.72 

282  +  79,12 

1285.22 

282  +  62.46 

1285.41 

282  +  46.79 

1285.43 

282  +  32.98 

1285.06 

282  +  20,12 

1284.80 

282  +  09.1  1 

1284,72 

281  +  98.11 

1284.43 

281  +  87,10 

1284.60 

281  +  76.09 

1283.30 

281  +  59.61 

1284.06 

281  +  43.10 

1284.64 

281  +  27.20 

1284.88 

280  +  95.40 

1283.96 

280  +  63.60 

1283.16 

280  +  31.80 

1283.37 

280  +  00.00 

1283.20 

tation 


289  +  03.59 


288  +  70.26 


288  +  36.92 


288  -H  03.69 


287  +  61.97 


287  +  20.35 


286  +  78.74 


286  +  37.12 


285  +  35.50 


285  +  53.88 


285  +  12.26 


284  +  70.64 


284  -f-  29.03 


283  +  87.41 


283  +  46.79 


283  +  12.46 


282  +  79.12 


282  +  62.46 


282  +  45.79 


282  +  38.08 


282  30.38 


282  +  25.25 


281  +  76.12 


281  +  75.50 


281  +  71.89 


281  +  65.46 


281  +  59.22 


281  +  50.18 


281  +  39.32 


281  +  23.85 


280  +  95.40 


280  +  63.60 


280  +  31.80 


280  +  00.00 


1293.15 


1292.56 


1291.70 


1291.26 


1291.14 


1290.76 


1290.12 


1 289.48 


1288.95 


1 288.52 


1288.09 


1287.63 


1287.11 


1286.59 


1286.U2 


11285.64 


1285.18 


1284.86 


1284.72 


1284.68 


1 284.48 


1283.94 


1291.17 


1281.01 


1280.98 


1282.37 


1283.48 


1284.00 


1284.04 


1284.36 


1283.99 


1283.92 


1284.11 


1 283,91 


1283.47 


tation 


291  +  54.94 


90  +  70.94 


90  -f  28.94 


90  +  07.94 


89  +  66.29 


89  +  45.64 


289  -I-  03.59 


288  +  70.26 


288  +  36.92 


288  +  03.59 


287  +  61.97 


287  +  20.35 


286  +  78.74 


286  +  37.12 


2tib  +  96.50 


285  +  53.88 


28b  +  12.26 


284  +  70.64 


284  -h  29.03 


283  +  87.41 


283  +  46,79 


283  +  12.46 


282  +  79,12 


282  +  62.46 


282  +  46.79 


282  +  32.98 


282  +  20.12 


282  +  09.11 


281  +  98.11 


281  +  87.10 


281  +  76.09 


281  +  69.61 


281  +  43.10 


281  +  27.20 


280  +  96.40 


280  +  63.60 


280  +  31.80 


280  +  00,00 


1298.48 


1297.06 


1296.21 


1296,72 


1296.72 


1294.41 


1293.29 


1292.65 


1292.43 


1293,67 


1293.63 


1293.39 


1292.70 


1291.94 


1291.22 


1290.90 


1290.62 


1290.29 


1289.96 


1289.54 


1288.98 


1288.33 


1287.68 


1286.43 


1284.79 


1284.66 


1284,60 


1284.37 


1284,29 


1284.16 


1284.62 


1284.19 


1284,30 


1285.41 


1285.40 


1284.24 


1283.29 


1283.71 


1283.91 


1283.64 


able  5  (Concluded) 


293  +  64.94 

1303.62 

293  +  22.94 

1302.61 

292  +  80.94 

1301.12 

292  +  38.94 

1301.30 

291  +  96.94 

1 300.46 

291  +  54.94 

1298.66 

291  +  12.94 

1297.17 

290  +  70.94 

1295.99 

290  +  28.94 

1296.10 

290  +  07. 

1294.94 

289  +  86.94 

1294.44 

289  +  66.29 

1293.10 

289  +  46.64 

1292.63 

289  +  03.69 

1292.45 

288  +  70.26 

1291.31 

288  +  36.92 

1290.36 

288  +  03.59 

1289.13 

287  +  61.97 

1283.48 

287  +  20,35 

1288.16 

286  +  78.74 

1287.91 

286  +  37.12 

1287.46 

286  +  96.60 

1206.78 

285  +  53.88 

1286.06 

286  +  12.26 

1285.40 

284  +  70.64 

1284.82 

284  +  29.03 

1284.31 

283  +  87.41 

1283.86 

283  +  45.79 

1283.41 

283  +  12.46 

1283.87 

282  +  79.12 

1284,44 

282  +  62.46 

1284,71 

282  +  46.79 

1284,91 

282  +  32.98 

1284.68 

282  +  20.12 

1284,33 

282  +  09,1 1 

1284,34 

281  +  98.1  1 

1284.19 

281  +  87,10 

1284.22 

281  +  76,09 

1283.97 

281  +  59.61 

1283.72 

281  +  43,10 

1284,09 

281  +  27.20 

1284.63 

280  +  95.40 

1283.57 

280  +  63.60 

1282.93 

280  f  31.80 

1282.61 

280  +  00.00 

1282.63 

tation 


289  +  66.29 


289  -»■  46.64 


289  r  03.69 


288  +  70.26 


288  +  36.92 


288  +  03.69 


287  +  61.97 


287  +  20.36 


286  +  78.74 


286  +  37.12 


236  +  95.60 


286  +  63.88 


285  +  12.26 


284  +  70.64 


284  +  29.03 


283  +  87.41 


283  +  46.79 


283  +  12.46 


282  +  79.12 


282  +  62.46 


282  +  46.79 


282  +  38.08 


282  +  30.38 


282  +  26.26 


282  +  20.12 


281  +  76.12 


281  +  75.60 


281  +  71.89 


281  +  66.46 


281  +  69.22 


281  +  60.18 


281  +  39.32 


281  +  23.86 


280  +  96.40 


280  +  63.60 


280  +  31.80 


280  +  00.00 


1293.28 


1292.82 


1292.54 


1292.16 


1291.66 


1291.27 


1290.92 


1290.45 


1289.85 


1289.26 


1288.73 


1288.29 


1287.84 


1287.36 


1286.81 


1286.29 


1286.76 


(1286.37 


1284.90 


1284.56 


1284.41 


1284.31 


1284.15 


1283.44 


1290.73 


1281.01 


1280.80 


1281.96 


1282.78 


1283.34 


1283.63 


1283.96 


1283.99 


1283.56 


1283.68 


1283.63 


1283.32 


tation 


292  +  38.94 


291  +  96.94 


291  s  54.94 


291  +  12.94 


290  +  70.94 


290  +  28.94 


290  +  07.94 


289  +  86.94 


289  +  66.29 


289  45.64 


289  -I-  03.59 


288  -h  70.26 


288  4  36.92 


288  -I-  03.59 


287  +  61.97 


287  +  20.36 


286  +  78.74 


286  +  37.12 


286  +  96.60 


286  +  53.88 


286  +  12.26 


284  +  70.64 


284  -I-  29.03 


283  +  87.41 


283  +  46.79 


283  +  12.46 


282  +  79.12 


282  +  62.46 


282  +  46.79 


282  +  32.98 


282  +  20.12 


282  +  09.1 1 


1301.34 


1 300.44 


1298.62 


1297.11 


1296.97 


1295.10 


1294.92 


1294.38 


11293.06 


1292.56 


1292.49 


1293.16 


1293.40 


1293.32 


1292.55 


1292.01 


12S1.47 


1290.99 


1290.67 


1290.16 


1289.76 


1289.34 


1288.90 


1238.37 


1287.77 


1286.66 


1284.87 


1284.64 


1284.41 


1284.31 


1284.12 


1284.03 


281  +  98.1 1 

1283.92 

281  +  87,10 

1284.18 

281  +  76.09 

1283.56 

281  4-  59,61 

1284.01 

281  +  43.10 

1284.88 

281  +  27.20 

1284.43 

280  +  96,40 

1283.31 

280  +  63,60 

1282.80 

280  +  31.80 

1283.33 

280  +  00.00 

1283.28 

NUMERICAL  MODEL 
COMPUTATIONAL  MESH 
DISCHARGE  19,000  CFS 
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Appendix  A 

The  Hydrodynamic  Modei: 
HiVEL2D 


HIVEL2D  is  designed  to  simulate  flow  typical  in  high-velocity  channels. 
The  model  is  a  finite  element  description  of  the  two-dimensional  shallow- 
water  equations  in  conservative  form.  The  model  does  not  include  Coriolis, 
boundary,  or  wind  effects  as  these  are  typically  not  important  in  high-velocity 
channels. 

Vertical  integration  of  the  equations  of  mass  and  momentum  conservation 
for  incompressible  flow  with  the  assumption  that  vertical  accelerations  are 
negligible  compared  to  horizontal  motions  and  the  acceleration  of  gravity 
results  in  the  governing  equations  commonly  referred  to  as  the  shallow-water 
equations.  The  dependent  variables  of  the  two-dimensional  fluid  motion  are 
defined  by  the  flow  depth  h,  the  x-direction  conqxinent  of  unit  discharge  />, 
and  the  y-direction  component  of  unit  discharge  q.  These  variables  are  func¬ 
tions  of  the  independent  variables  x  and  y.  the  two  space  directions,  and  time 
t.  Neglecting  free-surface  stresses  and  the  effects  of  Coriolis  force  as  these 
are  not  considered  important  in  high-velocity  channels,  the  shallow-water 
equations  in  conservative  form  are  given  as  (Abbott  1979)': 


.  if*  +  ifz 

dt  dx  dy 


*  H  =  0 


where 


Q  = 


h 

P 

Q 


(Al) 


(A2) 


'  References  cited  in  this  Appendix  are  included  in  the  References  at  the  end  of  the  main 
text. 
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p 


(A3) 


(A4) 


(A5) 


where 

p  =  uh,  u  being  the  depth-averaged  x-direction  continent  of 
velocity 

q  =  vh,  V  being  the  depth-averaged  y-direction  component  of 
velocity 

g  =  acceleration  due  to  gravity 

o„,axy,ay^,0yy  =  Roynoids  stresses  per  unit  mass  where  the  first  subscript 
indicates  the  direction  and  the  second  indicates  the  face  on 
which  the  stress  acts 

p  =  fluid  density 

Zc  =  channel  inveit  elevation 

n  =  Manning’s  roughness  coefficient 

Cg  -  dimensional  constant  (Cg  =  1  for  SI  units  and  Cg  =  1.486  for 
non-SI  units) 


A2 
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The  individual  terms  in  the  conservation  equations  are  as  follows: 

a.  Acceleration  force  per  unit  width: 


b.  Pressure  force  per  unit  width: 


c.  Body  forces  per  unit  area: 


d.  Bed  shear  stresses: 


Pg 


n^pyfp 


and  pg 


n^q^p^  +  qi 


O’ 


The  Reynolds  stresses  are  determined  using  the  Boussinesq  approach  of 
gradient'difhision: 


2k 


y» 


du 

'dx 

du  ^  dv 
dy  dx 


a  s  2v«— 

"  ‘dy 


(A6) 


where  v,  is  the  eddy  viscosity,  which  varies  spatially  and  is  solved  empirically 
as  a  function  of  lo^  flow  variables  (Rodi  1980). 
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A3 


V,  = 


(A7) 


where  Q  is  a  coefficient  that  varies  between  0.1  and  1.0. 

This  system  of  partial  differential  equations  is  solved  using  the  finite  ele¬ 
ment  method.  The  finite  element  approach  taken  is  a  Petrov-Galerkin  formu¬ 
lation  that  incorporates  a  combination  of  the  Galerkin  test  function  and  a 
non-Galerkin  component  to  control  oscillations  due  to  convection.  An  integra¬ 
tion  by  parts  procedure  is  used  to  develop  the  weak  form  of  the  equations. 

The  weak  form  facilitates  the  specification  of  boundary  conditions.  The  weak 
form  is  given  as; 


3jc  *  dy  ^ 
dO,  +  §4>i  [F/i,  +  F^n^) 


dQ 

dx 


(A8) 


where  the  variables  are  understood  to  be  discrete  values  and 
e  =  subscript  indicating  a  particular  element 
Q  domain 

=  0/  +  yj,  =  test  fimetion 
<t>i  =  Galerkin  part  of  the  test  function 
I  =  identity  matrix 

iPi  =  non-Galerkin  part  of  the  test  function 
riy)  =  n=  unit  vector  outward  normal  to  the  boundary  F, 


A 

B 


3F, 

ac 

3F, 

dQ 


(A9) 


A4 
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Natural  boundary  conditions  are  applied  to  the  sidewall  boundaries  through  the 
weak  statement.  The  sidewall  boundaries  are  "no-flux"  boundaries;  that  is, 
there  is  no  net  flux  of  mass  or  momentum  through  these  boundaries.  These 
boundary  conditions  are  enforced  through  the  line  integral  in  the  weak 
statement. 


Petrov-Galerkin  Test  Function 

The  Petrov-Galerkin  test  function  is  defined  <Berger  1993)  as: 
tA,  =  +  <fii 

where 

d(l>,  A  d0;  A 

where  /3  is  a  dissipation  coefficient  varying  in  value  from  0  to  O.S,  is  the 
linear  basis  function,  and  Ax  and  Ay  are  the  grid  intervajs.  A  detailed 
explanation  of  this  particular  test  function,  in  particular  A  and  B,  is  given  in 
Berger  (1993). 


Shock  Capturing 

Because  a  lower  value  of  /3  (/3  =  /S^)  is  more  precise,  a  large  value  of  0 
(0  =  =  0.5  where  0gnf  and  0^^  are  the  Petrov-Galerkin  parameters  for 

smooth  flow  and  for  shocks)  is  applied  only  in  regions  in  which  it  is  needed. 
HIVEL2D  employs  a  mechanism  that  detects  shocks  and  increases  0  automati¬ 
cally.  Therefore,  0gff  is  implemented  only  when  needed  as  determined  by 
evaluation  of  the  element  energy  deviation.  In  a  similar  manner,  the  eddy 
viscosity  coefficient  C  varies  from  Csn  to  Csh,  the  effect  being  that  eddy 
viscosity  is  increased  only  in  areas  of  greatest  element  energy  deviation. 


Temporal  Derivatives 

A  finite  difference  expression  is  used  for  the  temporal  derivatives.  The 
general  expression  for  the  temporal  derivative  of  a  variable,  Qj,  is: 

m*l  AT  f  • 

-a  g"'  *(!-«) 


(AlO) 

(All) 
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where y  is  the  nodal  location  and  m  is  the  time-stqp.  An  a  equal  to  1  results 
in  a  first-order  backward  difference  approximation,  and  an  a  equal  to  1.5 
results  in  a  second-order  backward  difference  approximation  of  the  tetiqioral 
derivative. 


Solution  of  the  Nonlinear  Equations 

The  system  of  nonlinear  equations  is  solved  using  the  Newton-Raphson 
iterative  method.  Let  be  a  vector  of  the  nonlinear  equations  computed 
using  a  particular  test  function  \pi  and  using  an  assumed  value  of  Qj.  Rj  is  the 
residual  error  for  a  particular  test  function  i.  Subsequently,  R,  is  forced 
toward  zero  as; 

A^;  -  -Rl  (A13) 

ae; 


where  the  derivatives  composing  the  Jacobian  are  determined  analytically. 
This  system  of  equations  is  solved  for  Aqj  and  then  an  improved  estimate  for 
is  obtained  from; 

(t;'  =•  e;  ^  h*  (aw) 


where  k  is  the  iteration  number.  This  procedure  is  continued  until  conver¬ 
gence  to  an  acceptable  residual  error  is  obtained. 


A6 
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The  San  Timoteo  Creek,  located  in  southern  California,  is  a  tributary  of  the  Santa  Ana  River  and 
drains  portions  of  the  San  Bernardino  and  Riverside  Counties.  The  existing  creek  has  the  capacity  to  protect 
the  surrounding  community  fiom  rqiproximately  a  20-year-frequency  flood.  The  proposed  chatmel 
improvements  will  provide  a  100-year  level  of  protection.  The  proposed  chaiuiel  design  within  the  reach 
studied  includes  a  sediment  basin,  a  concrete  weir  followed  by  a  converging  sidewall  chute,  a  compound 
horizontal  curve,  and  a  bridge  pier.  This  study  was  initiated  because  there  was  concern  as  to  the  adequacy  of 
a  one-dimensional  analysis  of  the  flow  conditions  within  the  channel  chute.  A  two-dimensional  analysis  was 
deemed  necessary  to  evaluate  the  chute’s  influence  on  the  flow  conditions  in  the  curve  and  the  curve’s  impact 
on  the  flow  conditions  at  the  bridge. 

The  two-dimensional,. depth-averaged  flow  model,  HIVEL2D,  was  used  to  simulate  the  flows  in  the 
high-velocity  channel.  This  model  was  chosen  because  of  its  ability  to  simulate  supercritical  flow  and  csq^ture 
shocks  such  as  oblique  standing  waves.  Simulation  results  indicated  that  the  proposed  San  Timoteo  Channel 
design  and  in  particular,  the  San  Timoteo  Canyon  Road  bridge,  will  convey  the  design  discharge  (100-year- 
frequency  event,  19,(X)0  cfs)  in  an  acceptable  manner. 
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